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Abstract 

We investigate the effect of boundary conditions on spin amplification in spin chains. We show 
that the boundaries play a crucial role for the dynamics: A single additional coupling between 
the first and last spins can macroscopically modify the physical behavior compared to the open 
chain, even in the limit of infinitely long chains. We show that this effect can be understood in 
terms of a "bifurcation" in Hilbert space that can give access to different parts of Hilbert space with 
macroscopically different physical properties of the basis functions, depending on the boundary 
conditions. On the technical side, we introduce semiclassical methods whose precision increase 
with increasing chain length and allow us to analytically demonstrate the effects of the boundaries 
in the thermodynamic limit. 
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I. INTRODUCTION 



Quantum state transfer through spin chains has attracted considerable attention starting 
with a seminal paper by Bose [1]. In such a scheme, an initial quantum state is prepared 
on a spin at one end of a chain, whereas all other spins are in a pre-defined state, say all 
pointing down. The system of coupled spins is then let to evolve freely, and after a certain 
time, long enough for a spin wave to propagate to the other end of the chain, the quantum 
state of the last spin is read-out. The last spin becomes in general entangled with the rest of 
the chain, and one therefore obtains a mixed state when ignoring the rest of the chain. Bose 
showed that the fidelity of such a transfer through an un-modulated spin chain with fixed 
nearest-neighbors Heisenberg couplings exceeds the maximum classically possible value for 
up to 80 spins. This work has been generalized in several directions [2-10]. Substantial 
effort was spent to increase the fidelity of the state transfer. Perfect state transfer was 
predicted for chains with couplings that increase like a square root as function of position 
along the chain towards the center of the chain, leading effectively to a rotation of a large 
collective spin [11]. Also, reducing the coupling between the terminating spins of the chain 
and the rest of the chain was shown to provide a recipe for perfect state transfer, at the cost 
of slowing down the transfer [12]. It was noted that arbitrarily high fidelity could also be 
achieved through dual-rail encoding [13] in two chains, even with randomly coupled chains 
[14]. 

Spin chains have also been studied in the context of spin amplification. Detecting single 
spins, and even more so, to measure their state is a formidable challenge [15]. Lee and 
Khitrin proposed a clever scheme of a "quantum domino", where an initially flipped spin 
leads to the propagation of a domain wall and ultimately the copying of the initial spin state 
onto a GHZ like state, (a| l)+/3\ t))| l)®^" 1 ' ->■ a\ l)® N +/3\ |) 07V . Kay noted a connection 
between quantum state transfer and spin amplification, which allowed to map insights from 
optimal state transfer to optimal spin amplification [16] and vice versa. Indeed, the same 
representation of the Hamilton operator in the two cases can be obtained by exchanging 
the couplings and basis functions at the same time. For spin amplification, one wants basis 
functions with a single domain wall and a hamiltonian which flips just the spin adjacent to 
the domain wall, inducing the domino effect. For quantum state transfer, the interesting 
basis functions all have a single excitation located on one of the N spins, and the hamiltonian 
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consists of nearest-neighbors exchange couplings. 

Recently, there has been interest in geometrical and topological effects in spin- networks. 
Quantum state transfer was extended to more complicated networks, in particular hy- 
percubes [11], and to quantum computing during the transfer [17]. It was shown that 
during the transfer along a chain, an arbitrary single qubit rotation can be performed by 
appropriately splitting and recombining the chain. Even two qubit gates can be performed 
by coupling incoming and outgoing chains that carry the qubits to and from a central chain. 
Also, topological quantum gates were proposed, in which the chain is closed to a ring, and 
threaded by an appropriate Aharonov-Bohm flux [17]. Very recently, topological effects 
were exploited for locally controlling the dynamics in spin-networks [18]. 

In this paper we study the influence of the boundaries on spin amplification in spin 
chains. One might think that the boundaries consisting of the two terminating spins should 
play a negligible role in the limit of very long chains, N — > oo. The effects of boundary 
conditions indeed vanish in the thermodynamical limit in most situations in physics. Well- 
known exceptions exist in the presence of long-range correlations, as for example exactly at 
a quantum phase transition [19]. But, surprisingly, it turns out that in spin chains, far away 
from any phase transition, the boundary conditions can drastically modify the dynamical 
behavior. The presence or absence of a single additional coupling between the last and 
the first spin can lead to macroscopically different time-dependent polarization even in the 
limit of arbitrarily long chains. We will demonstrate that a kind of "bifurcation" in Hilbert 
space can take place that explains the macroscopically different behavior. Depending on the 
boundary conditions, basis functions can be reached which may differ both in their physical 
properties, as well as in the scaling of the dimension of the basis with N. Furthermore, 
we will show that there are different ways of closing a chain to a ring, which not only 
drastically modify the physical behavior of the chain, but even lead to nonequivalent matrix 
representations of the hamiltonian, and different dynamics in the accessible part of Hilbert 
space. There is even a way of closing the chain such that the different topology is felt in one 
part of the Hilbert space, but not in another. 

To study the chains in the limit of very large N, we introduce an innovative semiclassical 
approach, the precision of which increases with increasing N. The method works well for 
dimensions of the relevant basis which scale linearly with N, and allows us to prove persistent 
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macroscopic differences in the physical behavior for iV — > oo. We start the analysis by 
reviewing the "quantum domino" system introduced in [20], and developing the semiclassical 
method at the example of linear chains with open boundary conditions (simply called "linear 
chains") in the following. 

II. LINEAR CHAIN 

A. Description of the system 

We consider a linear chain of Nl spins — ~ with nearest-neighbors interactions, whose 
hamiltonian is given by 

j n l -i 

H {L) = -^I t (U-i g> l fc+ i - <8> Z k+1 ) . (1) 

k=2 

Xk, Zk are Pauli operators acting on spin k. The coupling constant J will be set to J = 1/2 
throughout the paper, is an effective hamiltonian derived in [20] for a one-dimensional 

Ising chain of two level atoms with nearest-neighbors interactions, irradiated by a weak 
resonant transverse field. The hamiltonian has the physical meaning that when a spin Sk 
is surrounded by two spins (Sk-i and Sfc+i) of opposite sign, the operator Xk flips spin 
Sk- In the original model, this is achieved by the dependence of the resonance frequency 
of an atom on the state of its neighbors. If we consider the situation where the system is 
initially in the state |Ai) = | 4it2 • • • tN L ), i-e. the first spin is down, all the others up, 
the dynamic is restrained to evolve in a subset of the total Hilbert space of size Nl — 1. 
Initially couples |Ai) to | A2) = |4-i4-2T3 • • • Ta^)- Then in general the system couples 
|Afc) = |ii • • • ifctfc+i • • • tiO (k spins down, all the others up with k — 2, . . . , Nl — 2) to 
|Afe_i) and |Afc + i). Finally, at the end of the chain, |Ajv l -i) is reached, which is itself only 
coupled back to |Ayv L _2}- On the other hand, the state where all spins are initially up, is an 
eigenstate of of eigenvalue and is therefore stationary. Thus, a stimulated wave of 
flipped spins can be triggered by the flip of a single spin. In other words the system acts as a 
spin amplifier that amplifies the initial states \\) or |f) to a macroscopic polarization of the 
entire chain. Note that the restriction to the small subspace of A^ — 1 states |Ai) , . . . , IAjv^-i) 
is a consequence of the fact that the initially excited spin is at the beginning of the chain. 
A single excited spin in the middle of the chain leads to significantly different dynamics (see 
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section III). The states {|Afc)}, k G [1, . . . ,Nl — 1], form an orthonormal basis in which 
i s represented by 
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In the subspace considered, is therefore equivalent to a ID tight-binding hamiltonian 

1 N L -2 

with constant nearest-neighbors hopping, = — |A^} (Afc +1 | +h.c, whose eigenstates 



are Bloch waves, 
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The corresponding eigenvalues form a ID energy band, 

,p=l,...,N L -l. 
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The knowledge of the exact eigenvalues and eigenvectors of allows us to obtain an 

analytical expression of the propagator, 



Ui%{t) = (A k \e- 

N L -1 
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(5) 



In [20], this form of the propagator was used to study numerically the time dependent 
mean polarization. In spite of the exponential simplification of the problem in the subspace 
considered, compared to the dynamics in the full 2 Nl dimensional Hilbert space, each of the 
(Nl — l) 2 matrix elements of the propagator (5) still contains a sum of 0{Nl) terms. We 
now show that very precise approximations of can be obtained that only involve one 
or few terms. This allows us to obtain closed analytical expressions for the time dependent 



mean polarization. We propose two such approximations. The first one is based on an exact 
representation of in terms of Bessel functions. The second is of semiclassical nature. 
For a given t, both become the more precise the larger N^, i.e. the longer the chain. 



B. Representation of the propagator in terms of Bessel functions 

Considering that 4%(p,t) = F$ Q (-p,t), and that F*%(0 ,t) = F$ Q (N L ,t) = 0, we 
can double the summation range and evaluate the sum by Poisson summation, 

N L -1 

2 



O) = \ E F $,(p.*) 

p=-N L 
-j oo rN L -± 

= o E / e^F™(p,t)dp. ( 6 ) 



m =— oo " ^£ 2 



We have F^ Q (p,t) = F^ q (2Nl +p,t), i.e. all matrix elements have a 2A^ periodicity in 
p. Hence, the integrals are over a period of the integrand and we are allowed to shift the 
integration interval as we like. Setting x = np/Ni we arrive at the propagator in terms of 
the Bessel functions with the argument t, 

1 oo p2N L 

O) = 2 E / ^S(?. f )* 

m=— oo 

_ \ r / e 2iN L mx-it cos as |" e i(fc+fco)a; _|_ e -i(fc+fc )x"| ^ 
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Using J_ n (x) = (— 1)" J n (x) we can simplify the result, 

oo 

m=—oo 
oo 

-H) k ~ k0 E (-l) fco+Arim ^o +2 ^ m (t). (7) 

m=— oo 

For the case ko — 1 (i.e. the situation of a single initially flipped spin at the left edge that 
we are interested in), we get further simplification due to the identity J n _i(£) + J n+ i(t) = 
2(n/t)J n (t), 

U$(t)= I E (-lf Lm (k + 2N L m)J k+2NLm (t), (8) 

m=— oo 

for all 1 < < iVx — 1. Eq.(8) is an exact expression that satisfies the initial condition 

(0) = At times less or equal to Nl (a single propagation from left to right), 
can be well approximated by the term m = 0, 

U$(t)^ 2{ ~f l kJ k (t), (9) 

owing to the fact that \Jj~ {t)\ <C 1 for \t\ <C k and 3> 1. Fig.l shows an example of the 
time dependence of [T^i f° r = 20, & = 5. Visible disagreement of the numerically exact 
propagator and the single-Bessel function approximation develops only for t > 27. Each 
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FIG. 1: (Color online) Approximation with a single Bessel function [m = in Eq.(8), dashed red 
line) of U^{t) (continuous blue line) for k = 5, Nl = 20 

additional term in (8) increases the range of t where the Bessel approximation is valid by 
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Nl- E.g., with Nl = 20, k = 5, adding two more Bessel functions (terms m = ±1) we get 
the continuation of the plot in Fig.l, with a visible deviation of the approximation from 
the exact result only at t > 65 (see Fig. 2). Physically, these terms correspond to waves 
reflected from the right and left edges, which explains their irrelevance for times when the 
spin waves have not yet reached the corresponding boundaries. 

0.6 
0.4 
0.2 



-0.2 

FIG. 2: (Color online) Approximation with three Bessel functions (m = —1,0, 1 in Eq.(8), dashed 
red line) of uj^(t) (continuous blue line) for k = 5,Nl = 20 




C. Semiclassical propagator: WKB approximation 

In section II D we will attempt to obtain a closed analytical formula for the time dependent 
polarization. In order to do so, another simplification of [/£y is in order. In fact, the leading 



k-l 



term m — in (8), that gives U$(t) = 2( ~f - kJ k (t), valid for t < N L , lends itself to further 
approximation. From [21] (leading term of 9.3.15), we can obtain a "WKB approximation" 
for the Bessel function, 



Jk KB (t) = ^^==cos(0 fc (t)) for k < t, (10) 
(f> k (t) = fcarccos (j\ - Vt 2 - k 2 + ^ . (11) 

This approximation is commonly called the Debye approximation [22]. The name "WKB 
approximation" is motivated by the fact that the Bessel function is approximated by a sum of 
two exponential functions with slowly varying amplitude and phase, in close analogy to the 
well-known WKB approximation method. The corresponding approximation of Jjjf^ (t) can 



also be obtained from a semiclassical solution of the Schrodinger equation, i.e. a Van Vleck 
propagator. This will be presented in III C 3 for circular chains. The WKB approximation 
breaks down near the classical turning point, which corresponds here to k = t (see Fig.3). For 
k > t, 4>k{t) becomes complex and Eq.(lO) has to be replaced by an exponentially decaying 
function. In the vicinity of the turning point, a uniform approximation is called for which 
interpolates smoothly between the two regimes. Bessel functions can be approximated near 




FIG. 3: (Color online) Comparison of exact Bessel function (continuous blue line) and its "WKB 
approximation" (dashed red line), Eq.(lO) 



the turning point by an Airy function ([21]: 9.3.23) 



2^ 



(t-k) 



(12) 



This allows to determine precisely the domain in which the WKB approximation is valid, and 
where it is not. Fig. 3 shows the comparison for t = 100 between exact Bessel function and 
its WKB approximation. The two plots practically coincide up to the last maximum of </&(£) 
as function of k, observed at k m = t—A. After that the WKB approximation tends to infinity 
like (t — k)~* whereas the exact Jk{t) quietly goes to zero. The maximal allowed approach of 
k to t in the WKB estimate is given by the position of the last maximum of the Airy function 

Ai(— z) which takes place at z rn = 1.01879 and is equal to Ai(—z m ) = 0.535. The width of 

i 1 

the "bad" region where WKB is senseless is thus given by A = t — k m = (|) 3 z m pa 0.8t3. It 
grows with time, but slowly compared to t. 
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D. Mean polarization 



N L -1 



The total polarization of the chain is defined as P^ = Zj.. In a basis state |Afc), we 

k=l 

have the mean polarization (Afc| P^ |Afc) = Nl — 2k. If we put the system initially in the 
state | Ai), we obtain the time dependent mean polarization 



N L -l 



k=i 



(13) 



Now let us substitute the WKB approximation of the propagator into the expression of the 
mean polarization. We are allowed to do so only in the "good" interval where the summands 
are smooth functions of k and the sums can be approximated by integrals. Neglecting the 
exponentially small terms for k > t, we are thus led to 
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where [x] denotes the integer part of x. The first integral in (14) is proportional to t, 
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Let us make an estimate of the remaining terms. The second integral -4 J t *_ A • • • is of the 
order \/tA « t 2 / 3 which is small compared to t when t — > oo. The sum containing cos [2</>fc(i)] 
contains many summands of both signs which mutually almost cancel and can be neglected 
(rigorous estimates can be made using the Poisson summation formula). Finally, the last 
sum can be estimated as the number of summands A oc t 1 / 3 , times the maximal value of the 
summand (Nl — 2t) 2 (|) Ai(—z m ) oc t 1//3 . The result is again oc t 2 ^ 3 and can again be 
neglected compared to t when t is very large. Thus the estimates for the relative errors decay 




FIG. 4: (Color online) Comparison of mean polarizations: exact result (Eq.(13), continuous blue 
line) versus linear term of WKB approximation (Eq.(15), dashed red line), for Nl = 50 



like t" 1 / 3 as t — > oo, and Eq.(15) is therefore the leading behavior of the WKB approximation. 
While this approximation is rather crude, it nevertheless shows very good agreement with 
the exact result for times small enough not to allow spin waves to reach the end of the chain. 
E.g., for Nl = 50 and t — 40 the Bessel representation (which is basically exact here), gives 
P(t) = -17.922 whereas the WKB result is N L - jgt = -17.9601, i.e. reproduces the 
Bessel result with 3 digits. In Fig.4 we compare the exact mean polarization (13) with its 
semiclassical approximation (15) for a chain of 50 spins. From Eq.(15) we read off a constant 
rate of spin flips (corresponding here to the speed of propagation of a single spin-wave front) 
for times t < N L , \v {L) \ = jg. 



III. CIRCULAR CHAIN 



There are several ways of closing the linear chain (1) of Nl spins to a ring: 
1. By identification of SV L +i and S±: in this system S 2 , ■ ■ ■ , Sn l are subject to spin 

11 



flipping (see Eq.(16) below); 

2. By identification of Sjvl+i an d <Si an d introduction of an additional interaction term 
containing X\ in (1) that allows S\ to flip: in this system Si, . . . , Sn l are subject to 
spin flipping (see Eq.(19) below); 

3. By introducing a different coupling between the first spins and the last spins of the 
linear chain: in section III C, we will introduce a particular four-spin coupling between 
SjVi-i, Sn l , Si, S2 . 

We will now investigate the dynamics in these different chains. Si will designate always the 
spin chosen as the reference spin from which all the others are numbered (from S2 to Sn l 
turning clockwise on the chains). 

A. Closing by identification of spins Sn l +i and Si 

The simplest way of closing a linear chain consists in identifying the first and the last spin. 
When we talk about circular chains we will call the number of spins Nq- For the present 
subsection we consider Nq = Nl, such that Sat c +i = Si. This imposes a corresponding 
boundary condition on the wave function, but also implies an extra term in the hamiltonian 
since identifying Si with Sat c+ i allows to flip spin Sn c depending on the state of Sat c _i and 
Sn c +i = Si. Thus we have to add the term k = Nc to H[ L \ i.e. 



Note that Si still never flips as there is no term containing X\. In section III B we will 
consider the situation of complete rotational symmetry implemented already on the level 
of the hamiltonian, where the introduction of yet another flip term also allows Si to flip. 
The system described by Eq.(16) evolves inside a subspace whose Nq{Nc — l)/2 basis states 




(16) 



k=2 
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\Xi,j) (I € [1) — 1] and J < I) can be arranged in the form of a triangle, 

lxi,i) 

IX2,l) > 1X2,2) 

|X3,l) » 1X3,2) , |X3,3) 

|Xjv -i,i) ) ) \Xn c -i,n c -i ) • (17) 

We define \xi,j) as the state containing / consecutive spins down (all the others are up), with 
J — 1 spins down at the left of S\, Si down, and I — J spins down at the right of S\ (imagine 
St at the top of the chain). For instance, |xi,i) = \i1t2 ■ ■ ■ tiV c ), |X2,i) = II1I2T3 • • • tiVb), 

|X2,2) = |4lt2 • • • tN c -liN c ), |X3,l) = 14-142^4 •■• tiVc)) 1X3,2) = |4l42t3 • • • TiV c -l4iV c ) > 

|X3,3) = \i1t2 ■ ■ ■ 1~n c -2In c -iIn c ), etc. 

couples a state \xi,j) to at most four neighbors in the triangle: the one above, the 
one above to the left, the one below, and the one below to the right. Depending on the 

region in the triangle, not all four of the states exist. A state is coupled to exactly those 

(o) 

of the four states that do exist. In order to obtain a matrix representation of H\ , we 
order the states \xi,j) m the order |xi.i) , |X2,i) , IX2.2) , |X3,i) • • i- e - we introduce a single 
label L(J, J) = 1(1 — l)/2 + J- Therefore, top state 1 couples to states {2,3} , lower left 
hand corner state L = (Nc-tfiNc-z) _|_ j couples to state L — Nq + 2, lower right hand 
corner state L = ( N c- 1 K N c-^) _|_ j\j c _ ^ COU pl es to state L — Nc + 1, left border states L = 
+ 1 with I e [2, N c - 2] couple to states {L - I + 1, L + I, L + I + 1}, right border 
states L = 1 + I with / G [2, Nq — 2] couple to states {L — /, L + /, L + / + 1}, base 
states L = (JVc ~ 1) 2 (jVc ~ 2) + J with J = [2, N c - 2] couple to states {L — N c + 1,L — N c + 2}, 
and inside states L = /( - J ~ 1 ' > + J with / G [3, N c — 2] and 1 < J < I couple to states 
{L — I, L — I + 1, L + I, L + I + 1}. For example, the matrix representation of H{ for 
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N c = 5 (iMfcli = 10 basis states ) readg 

/oil 



o o o\ 
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1 1 



1 



1 1 



1 



1 1 




1 



1 1 



1 1 



(18) 



2 



1 



1 1 







1 







1 1 







1 1 




1 J 



The general structure of the matrix corresponding to a system of Nc spins can be easily 
derived. The matrix is real and symmetric, and we need to consider only the upper right 
triangle (H^ )l,l> with L' > L. To fill this part of the matrix we have to know for a state of 
the basis to which state below in the triangle (17) it is connected. Each state from the first 
line to the last but one of the triangle is only connected to two consecutive states below it 
in the triangle (the state directly below and the one below to the right). Therefore, in each 

(C) 

line of the upper right part of H{ there are only two consecutive non zero elements. Given 
two states \l and Xl+i °f a same line of the triangle, there is one common state among 
the states coupled to \l and the ones coupled to Xl+i- That explains that we observe 
staircase-structures extending over 1, 2, 3, . . . rows . Since there are Nc — 1 lines in (17) and 
each line with the exception of the last one gives rise to a staircase, there are altogether 
Nc — 2 staircases. From one staircase structure to the next, there is a shift of one line and 
one column, since the first state of a line is not connected to the last state of the preceding 
one. 

A brief remark is in order about the global structure of the distribution of non-zero 
elements over the matrix. Let us consider for each staircase-structure the first element, 
i.e. )i 2 , (#1^)2,4 > (-^1)4,7? • • • The coordinates (line, column) of these elements are 

given by the index L for the elements of the left border in (17), Xj = L = J(7 ~ 1 ' 1 + 1, and the 
element below it, Yj = Xj + 1 — J ^ J 2 I ~ 1 ^ + 1. If we plot the (Xj, Yj) points in a plan, we obtain 
a line whose slope X I+1 x = 1 + ~ converges for I — ^ 00 to 1 . This means that for Nc 3> 1 , 
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the matrix representation of H[ is very close to a matrix having four non zero diagonals 
parallel to the main diagonal, but a final curvature of the "off-diagonals" remains for all 
finite /. In the next section we see that straight lines of off-diagonal matrix elements over 
the entire matrix, parallel to the main diagonal, are obtained with the second way of closing 

(C) 

the chain. Lacking a viable analytical technique for diagonalizing H{ , we diagonalize the 
hamiltonian numerically, and derive the propagator. The results will be compared to those 
of the linear chain in section IV. 

B. Completely periodic chain 

Another way of closing the system is to impose periodic boundary conditions S^a+k — Sk, 
k G [1, Nc]- This amounts to subjecting also S\ to flip by adding an additional flip term X\ 
controlled by SV C and S2. The hamiltonian of this system thus reads 

, N C 

HP = -^X fc (l fc _i®l fc+ i-#*-i®3fe+i) • (19) 

k=l 

The physical important point is that in this system Si is flipped if its nearest-neighbors are 
of opposite signs (Sn c and S2), whereas in section III A, Si is never subjected to spin flipping 
and always remains in its initial down state. As a consequence, the number of basis states 
that are dynamically accessible is doubled. Indeed, allowing Si to be flipped, gives access to 
basis states related to the \xl) by flipping all spins. They form an ensemble of states having 
k consecutive spins down (k G [1, Ac — 1]), where in addition, for each k G [l,Nc — 1], 
we have the freedom of choosing the starting point of the sequence of flipped spins among 
the Nc spins. So the dimension of the accessible Hilbert space is now Nc{Nc — 1). It is 
interesting to note that the couplings between the basis states can be obtained by arranging 
the states on the vertices of (Nc — 1) nested A^-sided polygons as can be observed in Fig. 5 
for Ac = 5. A full line represents the coupling between two states. We label the states as 
can be seen in Fig.5: we follow a polygon until all of its states have been accounted for. 
Then we move to the next polygon further inside by moving to its first vertex on the center 
of the last link of the previous polygon, and so on. Using this numbering, we obtain a matrix 
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H. 



(c) 



V 



/ 



(20) 



(c) 

FIG. 5: (Color online) Coupling of states introduced by H\ , Eq.(19), for Nc = 5 

fen 

representation of if^ > with a simple structure. For example for Nc = 4 we obtain 

^0 1 1 0^ 
000001100000 
000000110000 
000010010000 
100100001100 
110000000110 
011000000011 
001100001001 
000010010000 
000011000000 
000001100000 
000000110000 
The general structure corresponding to a system of Nc spins can again be easily derived. 

Since the matrix is real and symmetric, we need to consider only the upper right triangle 

(c) 

(i?2 )l,u with V > L. We observe two neighboring off-diagonals, parallel to the main diag- 
onal, and a few non zero elements next to the (vanishing) main diagonal. The non zero terms 
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directly next to the main diagonal are the elements (H\ )kN c ,kN c +i with k E [1,-ATc — 2]. 
The first diagonal extends from {H^ )i,jv c +i 

to (fff>) 

Nc(N c -2),Nc(N c -i) (all elements equal 
to one), the second diagonal goes from (ifg )i,jv- c +2 to (H^ )nc(Nc-2)-i,n c (Nc-i) ( a U e l e ~ 
ments equal to one except the (kNc) th with k E [l,Nc — 3]). The fact that the two most 
straightforward ways of closing the chain lead to a relevant Hilbert space of dimension of 
O(Nq) instead of O(Nl), is rather interesting and will be explored further below. At the 
same time it would be interesting to find out whether there exist other ways of closing the 
chain that lead to dynamics resembling as closely as possible, the dynamics of the linear 
chain, and in particular to a relevant Hilbert space of dimension of 0(Nc)- In the following 
two subsections we present such a closure and profit from the machinery developed in section 
II to study the corresponding dynamics. 



C. Closing by a particular coupling between SjVc-Ij Sn C i Si, and S2 

We now study a particular way of closing the linear chain for which, as we are going to 
see, the system evolves in a Hilbert space whose dimension scales only linearly with the size 

(C) 

of the chain. We consider a system that has the same nearest-neighbors interactions as H% 
for spin 5*2 to spin SV c -i, and an additional particular coupling for spins Sjvc-i; £>n C i $1 
and S2 that leads to the following action: if spin Sn c -i and spin S 2 are the same, then spin 
Sn c and spin Si are flipped, otherwise nothing is done. The hamiltonian of such a system 
is 

JVo-l 

H 3 C = - ^ X k (lfc— 1 <8> lfe+i — Z k _i g) Z k+ i) 

k=2 

+^X Nc Xi (l Nc _il 2 + Z Nc _iZ 2 ) ■ (21) 

Note the small difference compared to H 2 ■ As before, we denote \xk) — 
|4i • • • sUtfe+i • • • tiVc) = the state where spins Si to S k are down and all the oth- 

ers up for k E [1 , N c - 1]. Then we define \x Nc +k-i) = |ti • • ■ Uik+i ■ ■ ■ In ) = | A fc>, 
k E [1 , No — 1], the state where spins Si to S k are up, all the others down. The Afc) are the 
mirror states of \A k ) in the sense that all spins are flipped. Let us suppose again that the 
system starts in the state \xi) = |4it2 • • • tiv c )- This state couples to \x2) = lii-Uts • • • tiV c ) 
and |x2(jv c -i)) = IT1T2 ••• tjVc-i4-JVc) • For general k E [2 , N c - 2], \xk) couples to \xk-i) 
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and |Xfc+i) an( i \XNc+k-i) = | A^) couples to A^_x) and |Afc+i). Thus, at the end of the 
chain, the system branches over into the mirror states, which in turn lead back to the orig- 
inal states and therefore give rise to a closed basis of 2(Nc — 1) states. We note, however, 
that neither the state with all spins up nor the state with all spins down is any longer an 

(C) 

eigenstate of H% , with consequences for the spin amplification discussed in Sec. IV. The 
matrix of the hamiltonian in the basis of the \xk), f° r k G [1, 2(Nc — 1)] is 

/ 1 1 \ 

10 1 "-. 
1 



H. 



(O 





1 



1 

1 1 
1 



(22) 



This is a well-known hamiltonian of a regular chain closed into a ring. Its eigenvalues and 
eigenvectors are given by 



Af> = cos 1; 



pn 



c-l 

1 



and 

2(JVo-l) 

IT E 



pirk 



e ^c-i \ Xk ) 



for p 6 [1, 2(Nc — 1)]. Finally, we obtain the expression of the propagator 



(23) 



U$ (t) 



(x k \e-^ C ' h \ Xk0 ) 



2(JV C -1) 

E 

P =i 

2(JV C -1) 

E 



( Xfc |$f)>e-^ C) ($f)| Xfc0 ) 



7 {c) 

k,ko 



(24) 



withF fc ( 2(p,t) 



i 



2(N C -1) 



e ^c- 1 e 



N c -1 
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1. Representation of the propagator in terms of Bessel functions 



Using that F^{p + 2{Nc — 1) ,t) = Fjfj^(p,t), we can use the Poisson summation 
formula to find 

00 p2(Nc-l) 

u $oW = E / ^ mp ^l( P ,t)dp 

m=— oo " 

00 

= E H) fc -* 0+am(JV °- 1) J»-* 0+am (JV -i)(t) . (25) 

m=— 00 

For the special case where the system starts in the state we obtain 

00 

= H)*- 1 E (-ir^Jk-i^miNc-Dit). (26) 



771= — OO 



5. Semiclassical propagator for small times ( t < N c — 1) 

The Bessel function vanishes exponentially fast if its index is large and exceeds its argu- 
ment. Therefore the infinite sum (26) reduces in fact to only a few summands significantly 
different from zero. In particular, if t < N c — 1 only a single term corresponding to m = 
or m = — 1 survives respectively for 1 < k < N c — 1 (states A) and N c < k < 2(N C — 1) 
(states A), 

U${t) ~ {-i^Jk-iit), 1 < k < N c - 1; 
U™(t) « (-l)^- 1 ^)*- 1 ^-*.!^), iV c < A; < 2(iV c - 1) . 

If 1 <§C A; (respectively 1 <C 2A' r c — fc — 1), these Bessel functions can be replaced by their 
WKB approximation (10). 



3. Van Vleck approach 

Before considering the polarization dynamics, we present an alternative and more 
physical approach for calculating the semiclassical propagator for k G [l,Nc — 1] and 
k < t + 1. The same problem as before can be approached by deriving the classi- 
cal hamiltonian associated to the system; the method and its generalization to chains 
with slowly changing parameters are described in the review [23]. The action of 
on a basis state, \xk) = 2 (IXfc-i) + can be rewritten in terms of shift 
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operators, if we upgrade the index Ho a continuous variable, i.e. \x(k)) = 

\ (\x(k — 1)) + \x(k + 1))) = \ (e~ a l dk \x(k)) + e 9 ^ 9k \x{k))). This is a reasonable approach 
in the limit of very long chains, Nq — > oo, which can be considered the classical limit. 
If we denote <fi = —id/dk the momentum canonically conjugate to the coordinate k, we 
can write \x(k)) = \ [e~ 1 ^ + e'*) \x{k)) = cos (0) \x{k)). In the classical limit we 

therefore obtain a corresponding classical hamiltonian H^ 1 = cos (</>). Since H% 1 does 
not depend on k, the momentum is an integral of motion connected with the energy by 
E = cte = cos(0). The trajectory of the motion follows from the canonical equation 
^(^) = dgty = — sin(<^)), hence k(t) = ko — tsin(0). The classical action is given by 

rk 

S(k,ko,t) = / <p{k)dk — E(k, ko,t)t = (k — ko)(f) — E(k, ko,t)t. We find two classical 

J kn 



>k 

paths, (pi(k,ko,t) = — arcsin(^^ fl ) and (p2(k,ko,t) = n + arcsh^^y^). Their associated 



classical energies are Ei(k 1 ko,t) = cos(0i(/c, k , t)) = j\/t 2 — (k — ko) 2 and E 2 (k,ko,t) 



cos(02(/c, k , t)) = —\\ft 2 — (k — ko) 2 . Finally, the associated classical actions are given by 
Si(k, k ,t) = (k-k )(/)i(k, k ,t)-E l (k, k ,t)t = -(k-k ) arcsin(^ a )-y / t 2 - (k - k ) 2 and 
S 2 (k,k ,t) = (k-k )Mk,k ,t)-E 2 (k,k ,t)t = (A;-A; )(7r + arcsin(^)) + v /t 2 - (k - k ) 2 . 
The semiclassical Van Vleck propagator is obtained by summing over all the classical paths, 



fSW = £ «p(' («.(*■*..«) - "•£)> • <»> 

a V 

where v a is the Morse index for the classical path a. For the case of the small times 
considered we have only two classical paths with V\ — — 1 and u 2 = and therefore, 

Ull(t) = Ul^\t) + ul^\t) (28) 



uYl Al \t) = J ; = exp(-zX) 

k ' ko W 2ix^/t 2 — (k — ko) 2 1 ^ 



with X = (k — ko) arcsin(^ a ) + \Jt 2 — (k — k ) 2 . It can be easily checked that (28) with 
ko = 1 is exactly the same as U^ B {t)^ = {-if' 1 J^f B {t) with J^J B (t) given by (10). 
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4- Mean polarization 

2(N C ~1) 

The operator for the total polarization is now = Z^. Expressed in the basis 

k=l 

of the {|Xfc)} ; the mean polarization is (xk\ \Xk) — {Nc — 2k) for k E [1, Nc — 1], and 
(Xfcl p(C) \Xk) = -(3N C - 2k - 2) for k e [N c , 2(N C - 1)]. If we put the system initially in 
the state we obtain the time dependent mean polarization, 

Nc-l 2 



(P (C) (t)) = (No -2k) U${t) 



k=l 
2(N C -1) 



- (3iVc-2A;-2) ujg(t) (29) 

k=N c 

As before, the WKB approximation of the mean polarization is obtained by replacing the 
exact propagator in Eq.(29) by its WKB approximation U^^ B {t)( c \ 

= J2(N c -2k)\uZ KB (t)^\ 2 



k<t+l 



(3N c -2k-2)\U^ KB (t)^\ 

fc>2(JV -l)-f+l 

J2(N c -2k)\j^ B (t)\ 2 



k<t+l 



£ (3N c -2k-2)\j%« o B _ 1) _ k+1 (t)\ 

fc>2(JV c -l)-t+l 



E ^cos^+lj 

~V 7T Jt 2 -(k-l) 2 



k<t+l 



s/e -(k- iy- 

_ £ l {3Nc 2k n) c«'[»P(Af C -l) -*+!.«]] 
fe >2(iv -i)-t+i n ^t 2 - (2(N C - 1) - k + l) 2 

Making the same kind of calculation as before and neglecting in the limit Nq ^> 1 the terms 
in the vicinity of the turning point (see the discussion about the estimation of remaining 
terms neglected for the expression of mean polarization in the case of the linear chain in 
section I.D), we finally obtain for the linear contribution of the mean polarization 

(P%Ut))^(Nc-l)-^. (30) 

We read off a rate of spin flips | = ^. It is reduced by 25% compared to the one obtained 
for the linear chain, |i/ L )| = There are now two propagating wave fronts superposed. 
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Each contributes with about 50% probability, but the additional front propagating to the 
left has at each time step a smaller number of flipped spins, which explains the reduced total 
spin-flip rate. We thus have the astonishing result that a single modified coupling at the end 
of the chain can drastically and macroscopically change the dynamics of the entire chain, 
even in the limit of arbitrarily long chains. This is a highly unusual situation, as in most 
physical systems boundary terms are negligible in the thermodynamic limit. Exceptions can 
exist for systems exactly at a phase transition [19], but this is clearly not the situation for 
our spin chains. 

In section IV we examine the differences in the dynamics for the three different ways 
of closing the chain in more detail and elucidate their physical origin. Before doing so, 
we would like to point out, however, yet another way of closing the chain, with the rather 
peculiar property of allowing a dynamics corresponding to different topologies of the chain 
depending on the subspace of Hilbert space considered. 

D. Existence of a circular chain dynamically equivalent to the linear chain 

Consider the situation where spins S2 to Sn c -2 are coupled through the same nearest- 
neighbor interactions as for H^ L \ but where we have an additional coupling for spins Sn c -2, 
Sat c -i, Sn c and S\ that leads to the following action: if spins Sn c _2 and Si are the same and 
different from Sn c , then spin Sn c -i is flipped, otherwise nothing is done. The hamiltonian 
of such a system is given by 

N c -2 

HP = - X k (U-i ® l*+i - Z k -x ® Z k+1 ) (31) 

k=2 

+ -Xjv _i (Itvc-2 ® li + Z Nc _ 2 <g> Z x ) (l Nc ® U - Z Nc <g> Z x ) . 

o 

In this system the states with all spins up or down are still stationary. One verifies that in 
the full Hilbert space spanned by the 2 N ° = 2 Nl computational basis states the matrix rep- 
resentations of and differ. Nevertheless, it turns out that the subspaces connected 
to the initial state |Ai) through the dynamics generated by and H± are identical for 
the two hamiltonians. Furthermore, the matrix representations of and in these 
parts of Hilbert space connected to | A x ) are identical, and thus lead to identical dynamics. 
We therefore have the interesting situation that the different topology of the circular chain 
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manifests itself only in a certain subspace of Hilbert space, whereas within the subspace 
relevant for the spin amplification problem one cannot distinguish the two hamiltonians 
through the dynamics which they generate, whatever the observable. 



IV. COMPARISON OF LINEAR CHAIN AND CIRCULAR CHAINS 

In the following subsections we compare the time evolution of the total polarizations as 
well as the total fidelities with respect to the two initial states of the first spin for the different 
hamiltonians presented above. We define the total fidelities F^ x \t) (X E {L, CI, C2, C3}) 
for H (L \ Hf\ Hf ] and H { °\ as 

F^\t) = F^ x \t)+F[ x \t) with (32) 

/ N X 

4 X \t) = (A (t)\ ( X) I tmXtr 
\m=l 
f N x 

F[ x \t) = (A 1 (t)\ I W(U I ) |A x (t)> , (33) 

\m=l 

with |A (0)> = |A ) = |ti ■■■tjvjt) and |Ai(0)) = |Ai) = \Uh •••tiv x )- These fidelities 
can be considered as amplification factors summed over the two initial basis states. They 
are a generalization of the usual fidelity considered in the spin transfer problem, where the 
fidelity of the state of the final spin with respect to the pure initial states of the first spin 

are considered [1]. As H^ L \ H[ C \ and H% conserve the basis state |A ) we will see that for 

(c) 

them the total fidelity is directly related to the average total polarization. However, 
does not conserve |A ). In this situation, the total fidelities are useful for comparing the 
overall performance of the different devices as spin amplifier. Note that all fidelities have 
initial value + 1 (JVc + 1), and are bounded from above by 2Nl (2Nc) for linear (circular) 
chains. 




A. Comparison of and H$ 

In Fig.6 we compare the time evolution of the mean polarizations for the linear chain 

(c) 

studied in section II with basis of size Nl — 1, and the circular chain H% studied in section 
III C with basis of size 2{Nc — 1)- We first consider the case Nq = Nl. We observe that 
the mean polarizations for the two types of chains oscillate with approximately the same 
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frequency. For t < Nl the mean polarizations decrease linearly, but with a difference in the 
slopes of 25% as predicted by Eq.(15) and Eq.(30). The question of the physical origin of 
the different behaviors arises. The simulations shown in Fig.6 were done for Nl = Nc = 51. 

Adding an additional local term to the hamiltonian for closing the chain would be ex- 
pected to lead at most to 1/Nl ~ 2% effect. Moreover, our analytical results show that the 
25% effect persists for Nl — > oo. A first guess about the origin of the different behavior 
might be the dimension of the involved basis sets. Since the dimension of the basis for 
({|Afc)} with k G [1, Nl — 1], seen in section II A) is half of the one of ({|A/-)} plus the 
corresponding mirror states { lAjt)} with k G [1, Nc — 1], as seen in section III C), we might 
think that this is a reason for the observed difference. In order to check this hypothesis, we 
recalculated the polarization for the circular chain for Nc = (A^ + l)/2 in which case the two 
bases have the same dimension. In order to compare the polarization of chains of different 
length, we subtract the initial polarization, i.e. we consider the change of polarization with 
respect to the initial state. As Fig.6 shows, we observe the same difference between the 
slopes of the two mean polarizations as function of time (for t < Nl/2) as before. Again, 
this is confirmed by the leading terms in Eq.(15) and Eq.(30), which are independent of Nl, 
Nc- This shows that the difference of basis sizes is not responsible for the different behavior 
of the two systems. The matrix representations of the two hamiltonians and i/3 , 

Eqs.(2) and (22), differ only by the two off-diagonal matrix elements Hi^ L _i = H^ L _i i, 
which equal 1 for the circular chain, but for the linear chain. Thus, one would expect at 
worst a correction of 0(Nl) to the eigenvalues and eigenstates of H% . This is confirmed 
by numerical diagonalization, which indicates that the average absolute difference between 
corresponding matrix elements of the two propagators decays even more rapidly, roughly 
as 1/A r |. Clearly, the slightly different matrix representations cannot explain the observed 
macroscopic differences either. Fig.6 also shows that the mean polarization of the circular 
chain oscillates approximately twice as fast for Nc = {Nl + l)/2 compared to Nc = Nl- 
This makes sense, as the oscillation period is set by the time it takes for the spin waves to 
reach the end of the chain. Also, for Nc = (iVx + l)/2, the amplitude of the signal is roughly 
half the one for Nc = Nl since the number of spins of the circular chain has been divided 
by 2 (up to one spin). 

All of this shows that the important difference in the two models lies in the physical 
structure of the basis functions. The single additional coupling in H% allows access to 
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a new part of Hilbert space (i.e. the additional basis states A*A). Since the |Afc) have 
opposite polarization compared to |Afc), their admixture reduces the total polarization signal 
compared to the linear chain. Thus, a kind of bifurcation takes place in Hilbert space, 
depending on the presence of the additional coupling, and the difference in the physical 
properties of the additionally accessible basis states leads to the macroscopically different 
behavior. 

In Fig.7 we show the total fidelities. For H^ L \ F^ L \t) behaves like the inverted po- 
larization (P( L \t)), scaled by a factor 1/2 and shifted by a constant. Indeed we have 
F( L \t) = F (L) (t) + if^f) with F (L) (t) = ^ + |(A (t)|P^|Ao(t)) and F$ L \t) = 
^ - i(Ai(t)|p( L )|Ai(t)). But since |A ) is stationary F (L) (t) reduces to Nl, and so 

F {L) (t) = 3 -N L - 1 -(P^(t)} (34) 

For H^ C) we observe that the total fidelity immediately deteriorates (whereas F^ L \t) in- 
creases initially), and remains much smaller than F^ L \t). This behavior results from Fq (t) 
which is not constant, as |Ao) is no longer a stationary state. The dynamics generated by 
i?3 starting from |Ao) evolves apparently in a subspace of dimension 2 Nc ~ 2 (verified numer- 
ically up to Nq = 9). In order to compare the total fidelities for the different chains on equal 
footing, all fidelities have been plotted in Fig.7 for chains of 8 spins. Not shown in Fig.7 is 
that if we wait a sufficiently long time, the total fidelity F^ C3 \t) can become large again for 
particular values of t. As an example, for Nq = 8 we obtain F^ C3 \290) ~ 13.09, close to 
the maximal possible value 16 for this chain. However, comparisons of the fidelity should 
be made for a given fixed time interval, as for sufficiently long times a spin-configuration 
arbitrarily close to a given one can be found [1]. 



B. Comparison of H^ L \ h[°\ and Hf> 

In Fig.8, Fig.9, and Fig. 10, we plot the behavior of mean polarizations (P {L) ), (P^) 
and (PP) for chains of different lengths. For Hp' 1 and > the mean polarizations were 
calculated by numerical diagonalization since generalizing the semiclassical approach to these 
hamiltonians would be rather cumbersome. For (P^), (pP) and (P 2 ) we observe that 
the slopes of the curves for sufficiently small times, where the mean polarizations behave 
linearly, are independent of the size of the chains. They equal \v^ L n — |^ ~ 1.69 (obtained 
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<P(t))-<P(0)> 




FIG. 6: (Color online) Comparison of mean polarizations (P(t)} — (P(0)} for linear chain (con- 

(C) (C) 
tinuous blue line), circular chain with Nc = Nl (dashed blue line), circular chain with 

(CD (O) 
Nc = (Nl + l)/2 (dotted blue line), and circular chains H\ (dashed red line), (continuous 

red line) with Nc = Nl = 51 




FIG. 7: (Color online) Comparison of total fidelities F(t) for linear chain (continuous blue 

(C\ (c) (C 1 ) 

line), circular chain (dashed blue line), and circular chains H\ (dashed red line), 

(continuous red line), with Nc = Nl = 8 



analytically, see above) in the case of H^ L \ 



numerically) in the case of H[ C ^ and H^ 3 . Therefore, the polarizations of circular chains 

(C) (C) 

H\ and H 2 reach their first minimum respectively about 2 and 1.28 times faster than 
the linear ones. Since |Ao) remains stationary for H[ C ^ and Fl\ , the total fidelities F^ 01 ^ 
and F^ C2 ^ follow the same law as Eq.(34) with corresponding mean polarizations (see Fig. 7). 
The origins of the different behaviors compared to are more complex here. First of 



r(C) 



3.39 and 



2.16 (obtained 
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<PL(t)>- <PL(0)> 




FIG. 8: (Color online) Mean polarizations (p( L \t)) - (P (L )(0)) for linear chains for different 

chain sizes (from N = 30 (top curve) to Nl = 130 (bottom curve) in steps of 20) 

(PCl(t))- <PC1(0)> 




FIG. 9: (Color online) Mean polarizations (P^ (t)) — (P^ (0)) for circular chains H\ for different 
chain sizes (from N = 30 (top curve) to Nl = 130 (bottom curve) in steps of 20) 

(C) (C) (C) 

all, the dimensions of the set of basis states connected by H\ , and H% already 

differ in their scaling with the length of the chains (scaling as 0(Nq) for H[ c \ HP and as 
O(Nc) for H^P). Secondly, the matrix representations are substantially different, and in 
fact nonequivalent, such that the spectra of the hamiltonians in the accessible Hilbert spaces 
are different. This can be observed in Fig. 11 where we compare the spectra for the different 
hamiltonians. The size of the different chains have been chosen such that the basis of H^ L \ 
#T\ H?> and H^> are of same dimension (7140 basis states with Nl = 7141, Nci = 120, 
Nc2 — 80, Nc3 = 3571, with Nci, Nq2, Nc3, the sizes of the circular chains corresponding to 
H[ c \ HP and HP). We can observe that the spectra for and roughly coincide, 
as do the spectra for HP and IlP . For and H^P we have seen in subsections II A and 
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<PC2(t)>-<PC2(0)> 




FIG. 10: (Color online) Mean polarizations (P 2 (t)) ~ (P2 (0)) f° r circular chains for 
different chain sizes (from N = 30 (top curve) to Nl = 130 (bottom curve) in steps of 20) 

III C that the energy levels are given by a ID tight-binding model with nearest-neighbors 
hopping, and thus follow the cosine laws in Eq.(4) and Eq.(23). For H[ C \ H2 we have in 
the limit of large Nq for most states essentially a 2D tight-binding model: All states inside 
the pyramid Eq.(17) are coupled to four nearest neighbors, and the same is true for the 
states on the vertices of all polygons in Fig. 5 for H% , with the exception of the innermost 
and the outermost polygons. However, the spectra differ substantially here from the usual 
cos(k x ) + cos(ky) spectra for a tight-binding model on a square lattice due to the unusual 

(C) 

geometry. For example, for Nq = 5, the lattice corresponding to H2 has the dihedral 
symmetry of a pentagon (C$ v , see Fig. 5). This symmetry leads to rather large degeneracies, 
in particular in the center of the spectrum, and to finite slopes of the dispersion relation at 
the edges of the spectrum, which are particularly visible for relatively small chain lengths 
(see Fig. 11). Finally, also the relevant basis states are once more different from those of 



V. CONCLUSION 

In this article, we have shown how the way of closing a linear chain into a circular chain 
can lead to significantly different dynamical behavior. A single additional coupling may 
open access to additional parts of Hilbert space which may differ in their dimensions (and 
even the scaling of the dimensions with the number of spins), the physical properties of 
the basis functions, or give rise to different, nonequivalent matrix representations of the 
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FIG. 11: (Color online) (a) Spectra for (blue squares, Nl = 7141), h[ (continuous red line 

with crosses, N c = 120), H^ C) (red circles, N c = 80), Hf ] (blue triangles, N c = 3571). The 
cosine dispersion law of Eq.(4) is plotted as continuous blue line. Only a subset of the ordered set 
of all eigenvalues E\ is plotted as function of I in order to show the dispersion relation, (b) Spectra 
for H( L \ H[ C \ H^ C) and Hf ] for N L = N C = 10. Same symbols as in (a). 

hamiltonian. This can modify the time-dependent polarizations on a macroscopical level, 

even though the hamiltonians only differ locally by one spin and its coupling to the rest of 

(O 

the chain. For example, the rates of change of polarization in the linear regime for Hf> 

and (see Eqs.(30) and (15)) differ by about 25%, even for arbitrarily long chains. In 

(c) 

this case, the number of basis states of is twice as large as the number of basis states 
of but remains of the same order O(Nl). When adjusted to the same dimension 

of the basis by changing the length of the chain, almost the same representation of the 
hamiltonian results, but the physical properties of the basis functions differ strongly. The 
rates of polarization change for and differ even from the one of H^ L \ by roughly 
a factor 2 and 1.28, respectively. Here, the number of basis states connected by H\ or 
scales as O(N^). The matrix representations of the hamiltonians have a substantially 
different structure, and also the physical properties of the basis states differ. It remains 
to be seen whether applications can be found which exploit the high sensitivity of the 
time dependent polarizations to the additional local couplings used to change the boundary 
conditions. 
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